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Measurement of the B c meson mass in the exclusive decay B c — > J /if) ir 
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(Dated: February 28, 2008) 

A fully reconstructed B c — > J/ip + tt signal is observed with the DO detector at the Fermilab 
Tevatron pp collider. Using 1.3 fb _1 of integrated luminosity, the signal is extracted with a sig- 
nificance more than five standard deviations above background. The measured B c meson mass is 
6300 ± 14 (stat) ± 5 (sys) MeV/c 2 . 



PACS numbers: 13.25.Hw, 13.20.He, 14.40.Nd, 14.40.Lb 
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B c mesons are predicted by the quark model to be 
members of the J p — 0~ pseudo-scalar ground-state mul- 
tiplet and to have zero isospin as the lowest-lying bound 
state of a bottom anti-quark and a charm quark [l[ . B c 
properties are of special interest because of this meson's 
unique status as a short-lifetime bound state of heavy 
but (unlike quarkonia) different flavor quarks. Measure- 
ments of its mass, production, and decay therefore allow 
for tests of theories under new approximation regimes or 
extended validity ranges beyond quarkonia. 

This analysis uses data collected by the DO detector 
between April 2002 and March 2006 at the Fermilab 
Tevatron pp collider operating at ^fs = 1.96 TeV. The 
data sample corresponds to approximately 1.3 fb _1 of 
integrated luminosity. At the Tevatron the most eas- 
ily identified decay modes of the B c have a J/ip meson 
in the final state and are either the semileptonic mode 
B c — > J/tplv (I = e, /i), a signal with much higher 
statistics and thus more suitable for lifetime measure- 
ments, or the hadronic mode B c — ► J/ifiir, more suitable 
for mass measurements given its fully exclusive recon- 
struction without the loss of an escaping neutrino. 

The CDF collaboration has published results on both 
decay modes @, Q , and has recently updated the B c mass 
measurement to M(B C ) = 6275.6 ± 2.9 (stat) ± 2.5 (sys) 
MeV/c 2 [H]. The present Letter is the first report by 
the DO collaboration of a fully reconstructed hadronic 
decay mode of this state. The measured lifetime [3j 
is consistent with the expectation of a shorter B c life- 
time than for other B mesons due to the presence of a 
charm quark. The B c mass has been predicted by var- 
ious theoretical models [f| and most recently [|| with a 
three-flavor (unqucnchcd) lattice QCD numerical algo- 
rithm that yielded the smallest theoretical uncertainty, 
with the result M{B C ) = 6304 ± 12 + o 8 MeV/c 2 , where 
the first error is the sum in quadrature of statistical and 
systematic uncertainties, and the second is due to heavy 
quark discretization effects. 

The DO detector is described elsewhere Q, and the el- 
ements most relevant to this analysis are the tracking de- 
tectors inside a 2 T superconducting solenoidal magnet 
and the muon detection chambers. For enhanced pre- 
selection efficiency, no specific trigger requirements are 
applied, but all events satisfy one of a suite of muon trig- 
gers, typically requiring at least one muon with trans- 
verse momentum {pr) above 3 GeV/c. The decay under 
study consists of a single detached secondary three-track 
vertex: B c — > J/-07T — ► p, + fi~n (charge conjugate modes, 
7r^" , are always implied) . Initial track selection extends to 
a pseudorapidity of \rj\ < 2.0 (where i] = — ln[tan(0/2)], 
and 9 is the polar angle with respect to the beam line), 
and rejects tracks with pt < 1.5 GeV/c. Selected final 
state tracks must satisfy quality requirements based on 
established minimal hit patterns and a goodness of track 
fit. Tracks identified as muons must have matching hits 



in all three layers of the muon detector. 

Event selection starts with the requirement of an 
opposite-charge muon pair that forms a common vertex 
and whose mass is consistent with that of the J /if) meson 
(between 2.85 and 3.35 GeV/c 2 ). There follows a search 
for a third track that, together with the muons, must 
form a common vertex with y 2 < 16-0 for the three de- 
grees of freedom. The J/ ip candidate must have px > 4 
GeV/c, and the third particle is assigned the pion mass. 
Thus formed, the B c meson candidate is required to have 
p T > 5 GeV/c. 

Further B c candidate selection places constraints on 
quantities that proved to be strong discriminators against 
combinatoric backgrounds. The impact parameter (IP) 
significance of any particle, reconstructed either from 
a single track or a combination of tracks, is 7 S j g = 
yJ[e T /<T{t T )] 2 + [ei/(T(e L )] 2 , where e T (e L ) is the trans- 
verse (longitudinal) projection (with respect to the beam 
direction) of the track IP relative to the pp interaction 
vertex, or primary vertex, and a is the associated un- 
certainty. The primary vertex is determined event by 
event using a method described in Ref. The trans- 
verse decay length significance of a decay (or secondary) 
vertex is S xy = L xy /a(L xy ) where L xy is the distance 
separating that vertex from the beam line. The pointing 
cosine, C xyi measures the alignment between L xy and 
the transverse momentum direction of the decaying can- 
didate particle. The isolation X of a B c candidate is 
defined as the ratio of two pt sums: that from the three 
candidate tracks, divided by that from all tracks with pt 
above 0.3 GeV/c whose momenta are lying within a cone 
of radius A1Z = J (A77) 2 + (Ac/)) 2 = 0.5, where A77 and 
A(f> are distances in pseudorapidity and azimuthal angle 
from the B c momentum axis, respectively. 

Throughout the background reduction process, a con- 
trol procedure is used that tests the effect of each dis- 
criminator against a well-understood signal sample, ei- 
ther reconstructed B ± — > J ji\)K ± candidates in data or 
candidates in a B^ — ► J/ipn^ simulated Monte Carlo 
sample. The latter isgenerated using EvtGen inter- 
faced with pythia |10j , followed by full modeling of the 
detector response with GEANT [ll[ and event reconstruc- 
tion exactly as in data. 

J/ifi candidates are mass constrained, i.e., their daugh- 
ter muon momenta are corrected to yield the PDG [12] 
mass value. When the third track is assumed to be a 
kaon, a clean, high-statistics B ± signal in invariant mass 
is observed in the data. This decay has a topology sim- 
ilar to the B c signal and is used as a reference in an 
initial round of selection cuts shown in Table U as Stage 
1. Here the S ± signal and sideband regions are used as 
efficiency and rejection indicators of where to set selec- 
tion thresholds. The S ± study region extends from 4.98 
to 5.58 GeV/c 2 in invariant mass, and the signal region 
is approximately ±2<r wide from 5.20 to 5.36 GeV/c 2 . 
Individual cuts are required to be about 95% efficient, 
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TABLE I: Discriminators and their values at the two selection 
stages (see text), (71") is introduced only for the second 
stage and represents the transverse momentum of the pion 
candidate with respect to the total B c candidate momentum. 



Discriminator 


Condition 


Stage 1 


Stage 2 




< 


3.5 


3.5 


Li g (n) 


> 


3.0 


3.5 


Sxy 


> 


3.0 


4.5 


Cxy 


> 


0.95 


0.95 


1 


> 


0.5 


0.64 


P t(tt) (GeV/c) 


> 


1.8 


2.2 


p^(n) (GeV/c) 


> 




1.5 


p r T cl M (GeV/c) 


< 




2.5 




M(uuji) [GeV/c 2 ] 



with typical background rejection of approximately 20%. 
The resulting thresholds are listed in Table HI 

However, there are differences between the B ± and B c . 
Due to the lower (by about 1 GeV/c 2 ) invariant mass and 
the longer (Mike versus c-like) lifetime of the B , back- 
ground reduction undergoes a second stage, in which the 
B c Monte Carlo is used to model the signal. This second 
selection stage (Stage 2 in Table UJ aims at re-optimizing, 
if needed, those cuts associated with B c specific decay 
properties. With the third track now assumed to be a 
pion, the range in invariant mass from 5.6 to 7.2 GeV/c 2 
is studied. A sub-range between 6.1 and 6.5 GeV/c 2 is 
treated as the B c signal search window, and its invari- 
ant mass distribution in data is kept blinded throughout 
the analysis. This sub-range is approximately ±3er (mass 
resolution as determined from simulation) wide, and cov- 
ers both the theory expectations for the B c mass @ as 
well as the observed values quoted in @, 0] . Data in mass 
sidebands outside this sub-range are used as a model for 
backgrounds and to quantify background rejection. Table 
U lists those selections that were re-optimized (or intro- 
duced, in the case of p^ 1 (tt)) in Stage 2, and summarizes 
their evolution between the two selection stages. At this 
stage there remain no dimuon vertices with more than 
one candidate for the third track, and no events with 
more than one B c candidate. 

From B c simulated events, the B c mass signal is found 
to be well-modeled by a Gaussian function with a width 
of 55 MeV/c 2 . The mass resolution of the B ± J/ipK* 
signal observed in the data under similar conditions, af- 
ter all selections have been applied, reproduces the same 
width when scaled by the ratio of the B^ and B c masses. 

The resulting J/-07T invariant mass is shown in Fig. Q] 
where a clear excess is seen near 6.3 GeV/c 2 . An un- 
binned maximum log-likelihood (UML) fit of the J/ipir 
invariant mass distribution is performed, where the sig- 
nal is modeled by a Gaussian function with width fixed to 
a value of 55 MeV/c 2 , and combinatoric backgrounds are 
modeled by a first-degree polynomial. The result of the 
UML fit is overlaid in Fig. [T] and yields a signal of 54 ± 12 
events and a B c mass value of 6300.7 ± 13.6 MeV/c 2 . To 



FIG. 1: J/ipir invariant mass distribution of B c candidates af- 
ter the final selection. A projection of the unbinned maximum 
likelihood fit to the distribution is shown overlaid. 

estimate the signal significance, the same fit is repeated 
under the assumption that no signal is present. From 
the negative log-likelihoods of the signal plus background 
and background-only hypotheses, the signal significance 
is extracted [13] as N a = {21n[£(s + b)/C{b)]} 1/2 = 5.2 
standard deviations above background. For another es- 
timate of signal significance, \ 2 fits to data (in the 40 
MeV/c 2 bins of Fig. Q} under both hypotheses produce 
an increase m fit x 2 of 27 units, again indicating N a = 5.2 
standard deviations above background. 

Possible biases and systematic uncertainties affecting 
the B c mass determination are estimated using both the 
B ± signal in the data and the B c signal in either the data 
or the simulation. Uncertainty assessments are made as 
these samples are re-fitted under various test hypothe- 
ses. Sources of systematic uncertainties are the event se- 
lection, the fitting procedure (input mass resolution and 
data modeling), and the reconstructed mass scale. 

The fitted mass values are examined in the simulated 
signal sample as the value of the pt(tt) threshold is var- 
ied from 1.9 to 2.5 GeV/c. No systematic mass bias is 
observed, but statistical fluctuations of ±4.0 MeV/c 2 are 
observed and assigned as a systematic uncertainty. Sim- 
ilarly, the p™ 1 lower threshold is varied between no cut 
and 2.0 GeV/c, and the resultant mass variation indi- 
cates a small upward mass bias of 0.5 MeV/c 2 for the 
cut value adopted with respect to the no cut case. The 
observed B c mass is corrected accordingly, and a 100% 
uncertainty is assigned to this correction. There is no 
indication of a bias in mass due to the upper p™ 1 limit. 

The values of the selection cuts that are not directly 
related to the kinematics of the third particle (the pion 
or kaon candidates in the B c or B ± cases, respectively) 
are varied within reasonable values. No mass biases are 
observed, and from the range of mass values obtained, a 
systematic uncertainty of ±2.5 MeV/c 2 is assigned due 
to the choice of these selection cuts. 
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TABLE II: Summary of systematic uncertainties in the B c 
mass measurement. 



Source 


Component 


Value (MeV/c 2 ) 


Selection 


7r kinematics 


4.0 




other 


2.5 


Data modeling 


mass resolution 


0.6 




background model 


0.5 




signal shape 


0.5 


Mass scale 




1.0 


Total 




4.9 



To assess the systematic uncertainty due to the uncer- 
tainty of the mass resolution, the width of the Gaussian 
is allowed to float in the fit. The width input is also 
changed from the nominal value of 55 MeV/c 2 to other 
fixed values in the ranee from 45 to 65 MeV/c 2 . From the 
variation of fitted mass results, a value of ±0.6 MeV/c 2 
is assigned to this uncertainty. 

The background model is changed from a first-degree 
polynomial to a second-degree and third-degree polyno- 
mial, and to an exponential function. From the result- 
ing change in mass observed, a systematic uncertainty 
of ±0.5 MeV/c 2 is assigned due to uncertainty in the 
background model. The signal model is changed from a 
single Gaussian to a double Gaussian function, and the 
resulting shift of 0.5 MeV/c 2 is assigned as a systematic 
uncertainty. 

Lastly, for an estimate of the mass scale uncertainty, a 
direct comparison is carried out between generated and 
reconstructed Monte Carlo masses, as well as between re- 
cent DO mass measurements of well-known B states and 



the world averages of their measurements 12] . From the 
observed range of mass differences, a systematic uncer- 
tainty of ±1.0 MeV/c 2 is assigned due to uncertainty in 
the DO mass scale for the B c decay. 

A summary of all systematic uncertainties in the B c 
mass measurement is shown in Table ILT1 The overall sys- 
tematic uncertainty is ±4.9 MeV/c 2 . The mass fit result 
of 6300.7 ± 13.6 MeV/c 2 is corrected by -0.5 MeV/c 2 
for the p^ 1 bias. The final result for the B c mass is 
6300 ± 14 (stat) ± 5 (sys) MeV/c 2 . 

In summary, using a dataset corresponding to 1.3 fb _1 , 
a signal for B c — > J/rfm has been observed with a signif- 
icance higher than five standard deviations above back- 
ground. The mass of the B c meson has been measured 
and found to be consistent with the latest and most pre- 
cise lattice QCD prediction Besides its relevance in 
the development and tuning of heavy-quark bound-state 
models, the B c sample described here, with added inte- 



grated luminosity, is expected to be used in the extraction 
of lifetime, relative branching ratio, and production rate. 
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